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FLIGHT INVESTIGATION AND THEORETICAL CALCULATIONS 
OF THE FUSELAGE DEFORMATIONS OF A SWEPT-WING 
BOMBER DURING PUSH-PULL MANEUVERS 
By Alton P. Mayo 

SUMMARY 


The results of fuselage -deflection measurements at three stations 
on the fuselage of a swept-wing bomber (Boeing B-47A) during 12 push-pull 
maneuvers are presented in the form of coefficients expressing the fuse- 
lage deflection due to the normal -acceleration loads , the pitching 
angular-acceleration loads, the pitching angular -velocity loads, and the 
combined effect of the zero-lift loads, fuselage dead weight, and landing- 
gear reactions. Comparisons are made between the experimental deflection 
coefficients and those obtained from theoretical calculations. An outline 
of the procedures used to obtain the experimental and the theoretical 
coefficients is presented. 


INTRODUCTION 


The design of a flexible airplane requires not only the determination 
of various aerodynamic and inertia load distributions that would occur 
under specified flight conditions but also the associated structural 
deformations . The computations are not only lengthy but are subject to 
numerous assumptions and the end result is sometimes in doubt. 

In order to secure an indication of the accuracy obtainable in each 
component of the calculations, as well as the accuracy of the end results, 
an extensive research program was undertaken in which a Boeing B-47A was 
used to obtain the pertinent experimental data. 

In one phase of the research program, deflections of various surfaces 
were measured for the purpose of comparing the measured deflections with 
the computed values. Some of the results of this phase have been presented 
in references 1 , 2 , and 3 , which give the deflection influence coefficients 
for the wing and an analysis of the wing deflections obtained in symmetri- 
cal flight and in rolling maneuvers. 
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The present report is concerned with the analysis of the fuselage 
deflections obtained under symmetrical flight conditions. These fuse- 
lage deflections were analyzed to obtain coefficients expressing the 
deflection due to the normal-acceleration loads, the pitching angular- 
acceleration loads, the pitching angular-velocity loads, and the combined 
effect of the zero-lift loads, fuselage dead weight, and landing -gear 
reactions. Comparisons are made between the experimentally determined 
deflections and those calculated theoretically. The procedures by which 
the flight data were analyzed for the various effects, and the methods 
used in the theoretical calculations are presented. 


SYMBOLS 

l section load, lb/in. 

n normal load factor measured at airplane center of gravity 

Z fuselage deflection measured from ground-zero position, positive 

downward, in. 

Z' fuselage deflection measured from no-load position, positive 

downward, in. 

q airplane pitching velocity, positive for nose pitching up, 

radians/ sec 

4 airplane pitching acceleration, positive for increasing positive 

velocity, radians/sec^ 

Sr? standard error of estimate of equation for Z, in. 

Subscripts : 


n 

derivative with respect 

to 

n 

for 

complete 

airplane 

4 

derivative with respect 

to 

4 

for 

complete 

airplane 

q. 

derivative with respect 

to 

q. 

for 

complete 

airplane 


0 zero lift 

1 combined landing-gear and dead-weight effects 
Matrix : 



column matrix 
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AIRPLANE AND TESTS 


The airplane used in the test was a Boeing B-V 7 A airplane. (See 
figs . 1 and 2 . ) The changes in the test airplane configuration from the 
standard airplane were the installation of an airspeed-measuring boom 
and fairing on the nose and an external canopy, housing the deflection- 
recording instruments, mounted atop the fuselage approximately at the 
intersection of the airplane center line and the wing 38 -percent-chord 
line . 


The flight-test data used in this paper pertain to push-pull maneu- 
vers flown during the B— 4-7A flight-research program conducted at the NACA 
High-Speed Flight Station at Edwards, Calif. 

The push-pull maneuvers were made at altitudes of 25,000 and 
30,000 feet with Mach numbers ranging from O .51 to 0.80. The maneuvers 
in the present report were made with airplane gross weights of approxi- 
mately 109,000 and 125,000 pounds and center-of -gravity positions of 
13 percent and 22 percent mean aerodynamic chord. The normal, load factor 
ranged from 0 . 3 g to 1 . 5 g and the pitching acceleration ranged from 0 .l 6 to 
-0.14 radian per second per second. The specific values of Mach number, 
altitude, airplane weight, center-of -gravity position, dynamic pressure, 
and the weight of fuel in each fuel tank are given in table I for each 
run analyzed in this report. These are some of the runs for which the 
wing deflections were analyzed in reference 2 . 


INSTRUMENTATION AND ACCURACY OF MEASUREMENT 


The instrumentation pertinent to this paper consisted of a pitch 
turnmeter, roll turnmeter, altimeter, airspeed indicator, normal accelerom- 
eter, aileron-position recorder, and an optigraph system for measuring the 
fuselage deflections. The roll turnmeter and the aileron-position recorder 
were used only to insure that the data selected were from push-pull maneu- 
vers with little or no roll. The accelerometer was located at 3^*25 per- 
cent mean aerodynamic chord, and the pitch turnmeter was located at 
27-19 percent mean aerodynamic chord. Corrections were made to the no rmal - 
accelerometer data for the small displacements of the instrument from the 
airplane center of gravity. All instruments were of the standard NACA 
photographically recording type with the exception of the optigraph system 
which was designed especially for the B-k^A airplane. 

The fuselage optigraph system consisted of three target lamps on 
the fuselage and an optical recording instrument located atop the fuse- 
lage approximately at the intersection of the 38 -percent -chord line of 
the wing and the center line of the fuselage. The recording instrument 
was a continuously recording camera mounted on a platform attached to 
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the wing carry-through section. (See fig. 3.) In order to facilitate 
the recording of the deflections with the optigraph in daylight, high- 
intensity infrared light sources were used in combination with infrared- 
sensitive recording film. 

The optigraph system was calibrated through the use of a calibration 
stick with several infrared lamps spaced at 6-inch intervals and held 
vertically at each target station. All inflight measurements were made 
with reference to the fuselage -drooped position with the airplane resting 
on its landing gear. 

The normal load factor and pitching angular -velocity values axe 
estimated to be accurate to ±0.01 and ±0.005 radian per second, respec- 
tively. The pitching angular-acceleration values were obtained from 
measurements of the slopes of the pitching-velocity trace and are esti- 
mated to be accurate to ±0.01 radian per second per second. The airplane 
weights listed in table I apply at the time of the maneuver and are esti- 
mated to be accurate to ±500 pounds . The fuel weights for each tank are 
estimated to be accurate to ±500 pounds per tank. The deflections are 
estimated to be accurate to ±0.15 -inch incremental and ±0.4 inch from 
ground-zero position. 


METHOD AND RESULTS 


In a symmetrical maneuver, the final fuselage-load distribution can 
be considered as resulting from the superposition of various types of 
load distributions. The fuselage structural deflections are assumed to 
be linearly related to these component loadings. In the analysis of the 
present paper, the final loading on the fuselage is considered to be 
composed of six components as follows : 

(1) A load distribution due to the fuselage angle of attack 

(2) A load distribution due to the wing -wash effects 

(3) A load distribution due to the loads induced by the wing 

(4) A load distribution due to the fuselage inertia loads (including 

horizontal and vertical tail) 

(5) A concentrated load at the c/4 of the horizontal tail due to 

the horizontal-tail aerodynamic loads 

(6) A concentrated load at the wing -attachment pins due to the wing 

aerodynamic and inertia loads 
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Analysis of the Flight Data 

In order to illustrate the type of data used in the analysis, typi- 
cal time histories of the target deflections and the associated airplane 
motions are shown in figure k. The variation of the fuselage deflection 
with normal load factor is shown in figure 5. 

The procedure, by which such measurements were reduced to coeffi- 
cients expressing the deflection of a point due to the various types of 
loads, makes use of the assumption that the small wing torsion effects 
which vary nonlinear ly with the wing loads may he considered negligible. 
Thus the loads on any station of the fuselage of the flexible airplane 
resulting from the inflight aerodynamic and inertia loadings is given by 
equation (AlO) of appendix A as 


1 = l o + z n n + z q4 + *q<l 


(1) 


Because the deflections at any point on the fuselage are linearly 
related to the loads, the corresponding equation for the deflections due 
to the inflight loadings is 


Z ’ = z o + V 1 + z q4 + z q<l 


( 2 ) 


The deflections as measured in flight were referenced to a ground- 
zero condition and therefore included the fuselage deflections due to 
fuselage dead weight, landing gear, and the wing dead weight. The meas- 
ured fuselage deflections may be expressed as 


Z = (z 0 + Zi) + Znii + Z q q + Z q q 


( 3 ) 


At this point it is sufficient to note that each Z coefficient in equa- 
tion ( 3 ) may involve all the six-component loads due to fuselage angle of 
attack, wing wash, tail -load effects, and so forth, enumerated previously. 
A more detailed interpretation of each Z coefficient is given in 
appendix A. 

In the data analysis the combined Z D and Zj_ deflections were 
obtained as one coefficient (Zo + Zq) which expresses the fuselage 
deflection change from the measuring reference position on the ground 
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to the zero-load -factor flight position in the air. These combined 
deflections are hereafter referred to as zero-lift-plus-droop deflections. 

Although it is relatively difficult to obtain the coefficients of 
equation (3) from flight data under conditions where n, q, or 1 are 
all varying, a least-squares method of data reduction was used in the 
present report to determine all of the coefficients simultaneously. 

During the tests, Mach number, dynamic pressure, weight, and center- 
of -gravity position were effectively constant for each run. Thus, for 
each target in each run, the deflections at the various stations may be 
represented in matrix notation 


jzj = (Z 0 + Zjjjlj + Znjn} + z q{<i}- + z q{l} 

where the columns j j- are corresponding values of Z, n, q, and q 
read from the flight records at 0.2-second time intervals during the run. 
The coefficients (Z 0 + Zi), Z n , Z<|, and Z^ were solved for by the 
method of least squares using approximately 25 data points per run. 
Least-squares procedure is described in detail in reference 3- 

The values of the (Z G + ZjJ, Z^, Z^, and Z^ coefficients obtained 

for each target in each run are presented in table II. Table III lists 
the coefficients, the standard errors of the coefficients, and the 
standard error of estimate S^, calculated by the method outlined in 
reference 3? for a typical run, that is, run 7 of flight 10. Typical 
fuselage -deflection curves formed by these coefficients are shown in 
figure 6. Plots of all the coefficients for each run are shown in fig- 
ure 7. In order to illustrate the ability of the coefficients to predict 
the fuselage deflections, a comparison of the measured fuselage deflec- 
tions with those determined by using equation (3) and the deflection coef- 
ficients is shown in figure 8 for flight 10, run 7* 


Theoretical Calculations 

In the calculations of the theoretical fuselage-deflection curves, 
the methods of reference 4 for the wing loads and wing wash were combined 
with the methods of references 5 and 6 for the loads induced on the 
fuselage by the wing and the loads due to fuselage angle of attack. The 
combination of these methods resulted in equations from which a simul- 
taneous solution of all the aerodynamic loads on the airplane were deter- 
mined. The procedure which was used is given in detail in appendix B. 
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The wing-structural stiffnesses used in the calculations were 
obtained from references 7 and 8. The fuselage structural stiffness 
was obtained from reference 9 with the exception of the fuselage sec- 
tion from station 742 to station 86l for which no data were available. 

The stiffness of this section was considered as an average of the sec- 
tion forward and rearward. The wing basic -weight distribution was 
obtained from reference 7 and was modified for the weight of the instru- 
ment installation. The fuselage -weight distribution was obtained from 
data on the XB-47 obtained from the manufacturer and was modified to 
include the weight of the instruments . 

The theoretical wing-twist matrices required for the calculations 
were obtained by the methods of reference 4. The wing -deflection matrices 
were obtained from references 1 and 2. The downwash angles along the 
fuselage center line were calculated by the methods of reference 4, 
appendix E. The wing section lift-curve slopes were determined from the 
data of reference 10 by the method of reference 4 and corrected for Mach 
number effects by the Prandtl-Glauert equation. 


Comparisons 

Comparisons are shown in figures 9, 10, and 11 between the experi- 
mental and theoretical deflections due to the normal acceleration loads 
and the pitching-acceleration loads for two runs, one at a fairly high 
dynamic pressure and weight and one at a lower dynamic pressure and 
weight. The comparisons were restricted to two runs because the theo- 
retical calculations were quite lengthy. Specifically, the comparisons 
made pertain to a Mach number of 0.60 at 25,000 feet, weight 
108,700 pounds, and a Mach number of O.76 at 50,000 feet and a weight 
of 125,400 pounds . 

The comparisons between the calculated and the experimental deflec- 
tions due to the normal acceleration loads are shown in figures 9 and 11. 
A complete solution, including all of the aerodynamic and inertia loads 
acting on the complete airplane and the theoretical deflections due to 
tail load and inertia effects only are shown. The deflections due to 
the tail load plus fuselage inertia were included in the figures to show 
that these effects are the major components of the fuselage deflections. 
The deflections due to the tail load plus fuselage inertia were obtained 
by using the theoretical tail load from the complete solution and the 
inertia loads , rearward of the front-spar attachment pin, were obtained 
from the fuselage -weight distribution. Details of the theoretical pro- 
cedures are given in Appendix B. The experimental data points shown in 
the figures are plotted incrementally from a line through the deflection 
curve at fuselage stations 425 and 515*5 which is where the wing bay, on 
which the optigraph was mounted, is attached to the fuselage. In fig- 
ure 9 ) the fuselage moment distribution is shown in order to illustrate 
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the shape of the distribution and the general magnitude of some of the 
components which combine to form the final moment distribution. 

Similar comparisons for the fuselage deflections due to pitching 
acceleration are shown in figures 10 and 11. Again the moment distri- 
butions are included for illustrative purposes. The experimental data 
points are again plotted incrementally from a line through the theoreti- 
cal curve at the wing -attachment stations. 


DISCUSSION 


Deflection Coefficients 

Independent of any theoretical comparisons, the validity of the 
deflection coefficients obtained from the least-squares reduction of the 
flight data was substantiated by the facts that the coefficients : 

(1) Reestablished accurately the flight measurements from the 
measured airplane motions . 

(2) Established fuselage -deflection curves of expected shapes and 
magnitude due to the flight variables. 

( 3 ) Showed erratic behaviors only when the deflections represented 
by the coefficient were approximately equal to the optigraph measuring 
accuracy . 

The agreement shown in figure 8 between the measured fuselage deflec- 
tions and the deflections calculated by using the deflection coefficients 
and the measured airplane motions indicates that equation ( 3 ) which was 
chosen to represent the data includes all of the important variables . 

The shape of the deflection curve, formed by the (Z 0 + Zi) coeffi- 
cients shown in figure 6, is believed to be due to the manner in which 
the fuselage bends between the landing gear and rear of the landing gear 
when the airplane is on the ground. Since this deflection curve involves 
several unknowns, it is difficult to establish its general shape, but it 
can be anticipated that the fuselage dead weight when combined with the 
landing-gear and wing dead-weight reactions would probably produce a 
fuselage -deflection curve with several inflection points. Examination 
of figure 7 in conjunction with table I shows that the (Z Q + Zj_) coeffi- 
cients decrease with decreasing Mach number and dynamic pressure as could 
be expected from zero-lift considerations. Isolation of the various 
effects which are combined in the (Z Q + Zi) coefficients was not possible 
because of the undetermined shape of the Zi deflection curve. 
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From theoretical considerations of the loads involved in the Z n 

coefficients as shown in equation (A9) of Appendix A , it can be deduced 
that the Z n coefficients vary with the weight of fuel in each tank, the 

dynamic pressure, and the Mach number. Several attempts were made to 
determine these independent variations but a sufficient amount and range 
of data were not available. However, examination of figure 7 in conjunc- 
tion with the flight conditions of table I does show that the airplane 
with the heavier fuel weight produced the larger additional fuselage 
deflections. The Zn deflections in every run were in the direction 
expected from considering only the inertia loads . 

The Z^ coefficients also vary with Mach number, dynamic pressure, 

and fuel weight distribution. Unsuccessful attempts were made to deter- 
mine the independent variations of the Z^ coefficients with these 

quantities, but sufficient data were not available. Examination of fig- 
ure 7 in conjunction with the flight conditions of table I shows that the 
deflections due to pitching acceleration are largely affected by the 
amount of fuel in the tanks. The deflection variation from this compo- 
nent is, however, in the opposite direction expected from inertia con- 
siderations as the aerodynamic tail load associated with pitching accel- 
eration is sufficiently large to overcome the inertia effects. 

The Zq coefficients are rather small and erratic and, as shown in 

table III, there are also fairly large standard errors for these coeffi- 
cients. The inflight deflections resulting from pitching velocity would 
be only about 0.1 inch as predicted from the coefficients and the meas- 
ured pitching velocity. Thus, the fuselage deflections due to pitching 
velocity are of about the same magnitude as the measuring accuracy of 
the optigraph system. 

The experimental deflection coefficients were obtained with fuel 
conditions ranging from approximately 50 -percent capacity to nearly full 
capacity. Within these weight limits, the Mach number and dynamic pres- 
sure at which the maneuvers were performed are sufficiently spread so that 
the incremental fuselage deflections could be closely predicted from the 
given coefficients for almost any longitudinal maneuver. 


Theoretical Calculations and Comparisons 

The flexibility of the fuselage between the wing carry-through bay 
attachment pins was considered in the theoretical fuselage -deflection 
calculations and this section of the fuselage had a high rate of curva- 
ture as is shown in figures 9 , 10, and 11. Even though this fuselage 
section was considered flexible, the wing carry-through bay on which the 
optigraph was mounted was considered rigid and experienced a rotation 
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and translation as the attachment pins moved with fuselage bending. The 
optigraph was, therefore, assumed to be measuring from a rigid reference 
plane through the wing -bay attachment pins . 

The effect of wing bending on the wing carry -through-bay deformations 
and thus on the optigraph measurements was also considered negligible. 
Calculations indicated that the maximum displacement of the carry-through 
bay at the airplane center line due to wing bending was approximately 
0.15 inch for these maneuvers. 

The assumption that the optigraph measured incremental values from 
the plane of the wing carry-through bay on which the optigraph (see 
fig. 3 ) was mounted is fairly accurate as evidenced by the agreement in 
figures 9* 10, and 11. In order to establish precisely the measurement 
reference line, the theoretical methods were set up to give deflection 
values at stations only 25 inches apart over the entire fuselage. The 
results were plotted to a scale 10 times as large as that of figures 
10, and 11 and the line was established and transferred to these figures. 

The theoretical methods predicted the deflections due to the normal- 
acceleration loadings very well at the two conditions selected for com- 
parisons; M = 0.60, weight = 108,700 pounds, and M = O. 76 , weight = 
125,400 pounds. From the bending -moment plot (fig. 9)> it is evident that 
the normal-acceleration deflections are mainly an effect of inertia. The 
agreement shown between the deflections due to the tail load plus inertia 
and the deflections including fuselage airloads indicates that fuselage 
airloads have little effect on fuselage deflections. As previously men- 
tioned, in calculating the deflections due to tail load plus inertia the 
effect of the inertia of the nose section was neglected, thus making the 
moment distribution zero at fuselage station 425 as shown in the figure. 

The experimental and theoretical deflections due to pitching accel- 
eration also showed very good agreement. In the case of pitching accel- 
eration, as in the case of normal -acceleration effects, the deflections 
are due mainly to tail load plus inertia. In the case of pitching 
acceleration, however, the effect of the tail load exceeded the effect 
of inertia whereas in the case of the normal-acceleration deflections, 
the effect of tail load was small as compared to the effect of inertia. 

No attempts were made to calculate the (Z 0 + Zj_) coefficients due 
to the difficulty in determining the zero-lift loads and the ground- 
reaction forces at the various landing gears, nor were calculations made 
for the deflections due to pitching velocity because of the smallness of 
the measured deflections and the somewhat erratic behavior as indicated 
by the large standard errors. 

In addition to the good agreements between experimental and theo- 
retical deflections, the wing-root angle of attack, the tail load, and 
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the wing loads which were obtained directly from the theoretical solution 
also showed good agreement with the experimental data obtained from ref- 
erence 11 and from preliminary results of strain-gage measurements on the 
wing and tail of the test airplane as obtained in the Flight Research 
Division of the Langley Aeronautical Laboratory as shown in table IV. 

The slight differences between the calculated center of gravity and the 
experimental value pertaining to the maneuver may be partly due to the 
fact that in the calculations the center of gravity of each section of 
the wing was assumed to be at the 38-percent-chord line. 

The good agreement shown in the comparisons indicates that the fuse- 
lage deflections on this airplane may be calculated with a good degree of 
accuracy by existing methods from the known geometric structures and 
loadings . 


CONCLUDING REMARKS 


The inflight deflections of the fuselage have been analyzed by the 
method of least squares to obtain the component deflections due to the 
normal-acceleration loads, the pitching-acceleration loads, the pitching- 
velocity loads, and a combined loading due to the zero-lift loads, the 
fuselage dead-weight and landing-gear reactions. This type of analysis 
was validated by the consistency of the values obtained for the various 
coefficients and the ability of the coefficients to predict the fuselage 
deflections from the measured airplane motions . 

The experimental data were obtained with fuel conditions ranging from 
approximately 50 -percent capacity to full capacity. Within these weight 
limits the range in Mach number and dynamic pressure at which the maneu- 
vers were performed is sufficiently large that the incremental fuselage 
deflections could be closely predicted for almost any longitudinal 
maneuver . 

The comparison between the experimental fuselage deflections and 
those calculated with the fuselage weight and stiffness distributions by 
available aeroelastic methods indicates that adequate methods exist or can 
be devised for calculating the fuselage deflections on this airplane in 
quasi-static maneuvers. 

The fuselage deflections on the airplane are mainly due to the 
effects of tail load and fuselage inertia. The inertia effect was the 
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major component of the deflections due to normal acceleration, whereas 
the tail-load effect was the larger component of the pitching-acceleration 
deflections . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 16, 1956. 
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APPENDIX A 

FUSELAGE DEFLECTION COEFFICIENTS 


The fuselage loadings which are associated with the fuselage deflec- 
tion coefficients are described in detail in this appendix. This detailed 
knowledge of the loads involved is also necessary in the calculation of 
the theoretical deflection coefficients for purposes of comparison. 


Symbols Pertinent To Appendix A 
m airplane pitching moment, lb-in. 

L total load, lb 

l section load, lb/in. 

n normal load factor measured at airplane center of gravity 

xt longitudinal distance from c /4 of horizontal tail to airplane 

center of gravity, in. 

Z' fuselage deflection due to inflight loadings, positive down- 

ward, in. 

q. airplane pitching velocity, positive for nose pitching up, 

radians/sec 

q. airplane pitching acceleration, positive for increasing posi- 

tive velocity, radians/sec^ 

W weight of the airplane, lb 

a angle of attack, radians 

Oj. wing-root incremental angle of attack (measured from the 

n = 0 position), radians 

Differential operators : 

variation with wing-root angle of attack for flexible airplane 

aa r,i=0 with no inertia forces acting 


NACA m L56LO5 


variation of wing-root angle of attack with normal load factor 
for complete airplane with aerodynamic and inertia forces 
acting (q = q = 0) 

variation of load with normal load factor for flexible airplane 
with no inertia forces acting 

variation with normal load factor due to inertia loads plus 
aerodynamic loads caused by inertia deformations 

variation with pitching acceleration due to q inertia loads 
and aerodynamic loads caused by inertia deformations 

variation with pitching velocity due 
caused by pitching velocity (pitch 
station where a^. is defined) 

Subscripts : 

n differential d/dn for the complete 

4 differential d/dq for the complete 

q differential d/dq for the complete 

t pertaining to the tail 

fj loads on the fuselage 

f fuselage 

0 zero lift 

1 inertia effects 

w wing 

a additional loads 

Fuselage Load and Deflection Equations 

The loads on the fuselage which cause the deflections are divided 
into the following main parts: 


to aerodynamic loads 
axis for this case is 

airplane 

airplane 

airplane 


Ik 


daj. 

dn 


dl 

^=0 


d 

d-j_n 



_d_ 

dq 
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(1) The load imposed on the fuselage by the horizontal -tail aero- 
dynamic loads. This load is considered concentrated at the c/4 of 
the horizontal tail. 

(2) The net load imposed on the fuselage by the wing, which is a 
result of the wing aerodynamic and inertia loads. This load includes 
the effects of wing wash, induced loads from the wing and the wing- 
attachment -pin loads. 

( 3 ) The load on the fuselage due to fuselage angle of attack com- 
bined with the fuselage loads due to its own inertia. The fuselage- 
inertia loads are considered to include the inertia load on both the 
horizontal and vertical tail. 

These loads are a grouping of the fuselage loads enumerated in the 
section entitled Method and Results." This loading breakdown may be 
expressed in equation form as 


fuselage fuselage due to tail + ^fuselage due to wing + 


-fuselage due to fuselage 


Each one of these three main loads is treated as a result of 

(1) The loads on the flexible airplane with no inertia forces acting. 

(2) The inertia loads combined with the aerodynamic loads resulting 
from the inertia deformations. 

If the small variation in the center of pressure of the horizontal 
tail due to elevator deflection is neglected, then the aerodynamic tail 
load may be expressed by: 


L t = 


"V+fjO + 


[_^w±f_ 

\ da r,i=0 > 


K 


f^W 

l d,n I 


n + 


^w+f ) 

V d i4/ 


1 + 



-*- x t 


(ai) 


equation 


These moments are about the airplane center of gravity. In this 

drcy+f 

cU n 


is due solely to the aerodynamic loads arising from the 


inertia deformations since by definition of the center of gravity, the 
airplane inertia moment about the center of gravity due to load factor 
is equal to zero. 
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Thus the tail load may he expressed as 


L t “ 


L t,o + 


dl't, 


idcLj*^ 1=0 


°r 



dci j 


(A2) 


It should be noted that ^^4 is the aerodynamic tail load which 

diq 

balances the total airplane inertia moment and the aerodynamic effects 
resulting from the pitching-inertia deformations of the wing. 

The load on the fuselage due to the fuselage -aerodynamic and 
inertia loads may be expressed by 


dZ, 


'dZ-A /'dZ 


f > f ■ if ’° * T + v* 


dlu 


-)n + ( ( A 3) 


Also the loads on the wing are given by 


^w “ ^w.o 


AY • (§)- * fe> * (?) 


(A4) 


Since the loads on the fuselage due to the wing are linearly related to 
the wing loads, it follows that 


, ( di f,w \ 

f,« - *f,w , 0 + ferid'v ^ d i n / n 


/dZ 


di4; 


q. + 


( dZ 


f,w 


\ dq 


(A5) 


By considering the combined effects of equations (A2), (A3), and (A5), 
it can be seen that the load at any station on the fuselage may be 
expressed by 


Z 


Z 0 + 


dZ 


\ da r,i=0 > 


K- + 





(A6) 


At all times during a maneuver the additional aerodynamic load on 
the airplane must balance the normal inertia forces plus any normal forces 
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resulting from the pitching-acceleration loads and the pitching-velocity 
loads so that 

0^ = nW + + Lg (A 7 ) 



where L 0 is the normal resultant of all of the additional aerodynamic 

ci 

loads on the airplane. Dividing equation (A7) by W yields: 


fe) 


“r 


— n + -rr + 


L q L q 


¥ T 


where S£L_ is the variation of n with a r for the airplane with q 

dtXj, 

and q = 0 . Thus 


a,. 


= L + is + 

l W W /\dn / 


(A8) 


Substituting equation (A8) into (a 6) yields 


Z — Ir, + 


dZ 


Idn-i 


1 = 0 / 


r dl \ 

V d l n / 


n + 


dZ V L 4 l 


dZ 


, dn i=0/ \ w J \ d lV 


<1 + 





dZ 


dZ 
1 dq/ 


(A9) 


Each term of equation (A 9 ) involves all of the component loads 
enumerated previously. 


dZ 


^ 1=0 


The term 


is defined as 
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Thus the term — is the variation with normal load factor of 
<^ 1=0 

the additional fuselage load due to the flexible airplane with no inertia 
forces acting. 

The loads on the fuselage are therefore expressible by 

Z = Z Q + Z n n 4- Z^q + (AlO) 

and since the deflections are linearly related to the loads : 


Z T - Z Q + Z n n + Z^q + Z^q 


(All) 
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APPENDIX B 

OUTLINE OF THE METHODS USED TO OBTAIN THE 
THEORETICAL DEFLECTION COEFFICIENTS 


In the calculation of the theoretical deflections of the fuselage, 
the effects of the fuselage inertia loads, the fuselage aerodynamic 
loads, the tail load, and the forces at the wing -attachment pins must 
be considered. The final fuselage aerodynamic -load distribution is 

considered to he the sum of the components due to (l) the fuselage angle 
of attack, (2) the wing wash, and (3) the wing induction effects. Thus, 
it is evident that the wing, fuselage, and tail loads are interdependent; 
therefore, the solution for the fuselage loads due to normal acceleration 
must be obtained from an equilibrium equation which involves all of the 
airplane normal -acceleration loads . A similar procedure is required to 
determine the fuselage loads due to pitching acceleration. 

The solution for the loads in the present case was obtained by a 
method somewhat similar to that of reference A which involves solving 
the wing load equations simultaneously with the total airplane force and 
moment equations which are expressed in terms of the wing loads, the 
fuselage loads, and the tail load. This solution yields directly the 
wing loads, the wing-root angle of attack, and the tail load. Using 
these knowns, all of the fuselage loads can be obtained. 


Symbols Pertaining to Appendix B 
A fuselage cross-section area 

a maximum radius of fuselage 

b wing span 

b e effective wing span for calculating wing-induced drag 

(b e = 0.9 X b from ref. 12) 

c wing section streamwise chord measured at center of the section 

c wing mean aerodynamic chord 

c d,o wing section drag coefficient at cj = 0 

c^ wing section lift coefficient 
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~ dCm / da 


for inboard nacelle 


/m a,e 


dC m /da for outboard nacelle 


El 

F 

h 

H 

l t 

H 

Lf 

M 

m 

niQ 

q* 

S 

T 1 

T 2 

x 

y 


section pitching -moment coefficient 
fuselage structural stiffness, lb-in. ^ 
concentrated load, lb 

semispan of bound leg of horseshoe vortex 

height of plane of wing above fuselage center line 

distance from c/k of horizontal tail to airplane center of 
gravity 

tail load, lb 

load on fuselage, lb 

bending moment , lb -in . 

airplane pitching moment, lb -in. 

wing section lift-curve slope, per radian 

dynamic press lore, lb/sq in. 

wing area, sq in. 

thrust of inboard nacelle, lb 

thrust of outboard nacelle, lb 

longitudinal ordinate 

lateral ordinate (measured from and perpendicular to fuselage 
center line) 

angle of attack of fuselage, radians 


a i,l 


streamwise angle of twist at inboard nacelle due to wing 
inertia, radians 


°±,2 


streamwise angle of twist at outboard nacelle due to wing 
inertia, radians 
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wing -root angle of 

attack, 

radians 


“ 1,4 

wing twist due to 

r— 1 

II 

• CJ* 

inertia loads, 

radians 

(X l,n 

wing twist due to 

1 — 1 

II 

inertia loads. 

radians 

a 

angle of attack, radians 




81 1 deflection of inboard nacelle due to inertia 

5 j_ 2 deflection of outboard nacelle due to inertia 


Matrix symbols : 

[a] structural matrix expressing wing angle of twist caused by 

nacelle pitching moment (see ref. 4 ) 



deflections at center of wing sections due to combined wing 
aerodynamic and inertia loads at flight condition, in. 

matrix of wing section drag, lb 

matrix of lengths of bound segment of wing horseshoe vortices 



influence coefficient matrix expressing the wing deflection at 
inboard nacelle due to wing 0.25 chord airloads 



influence coefficient matrix expressing wing deflection at 
outboard nacelle due to wing 0.25 chord airloads 



matrix expressing induced loads on fuselage in terms of wing 
loads 



matrix expressing section loadings on fuselage in terms of 
fuselage section angle of attack 



structural matrix expressing deflection of selected points on 
wing in terms of quarter-chord airloads 



structural matrix expressing wing twist at inboard nacelle due 
to wing quarter -chord airloads 
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structural matrix expressing wing twist at outboard nacelle 
due to wing quarter-chord airloads 

wing-section lift (including all flexibility effects), lb/in. 
load on wing panels of span 2h, lb 


[ s i + s i] 


wing aerodynamic downwash matrix calculated from geometrical 
structure of wing and fuselage vortices 


[ S 1,d| 

downwash matrix relating downwash at fuselage center line to 
wing loads 

M 

structural matrix expressing the wing streamwise twists 
wing quarter-chord airloads 

due to 

ra 

fuselage overvelocity matrix expressing angle of attack 
due to flow around fuselage (see ref. 4.) 

on wing 


matrix of lengths of fuselage sections 


H 

moment arms to airplane center of gravity, in. 


M 

deflection at selected points on wing due to wing and nacelle 
inertia 

[-] 

diagonal matrix in which diagonal elements are equal to 

1 

M 

wing section streamwise angle of attack in the flight condition 
as defined by equations (B20) and (B24) 

Matrix 

notation: 


{ } 

column matrix 


[ ] 

square matrix 


[°i 

diagonal matrix 


L J 

row matrix 
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Subscripts : 

wing quarter -chord airload 
iy. wing inertia 

FS wing front-spar attachment pin 

RS wing rear-spar attachment pin 

nac nacelles 

w wing 

T thrust of nacelles 

t tail 

D wing drag 

L c chordwise component of wing section lift vector 

f fuselage 


Fuselage Load Equations 

Since the aerodynamic loads on the airplane in symmetrical maneuvers 
are symmetrical about the airplane center line, only one -half of the air- 
plane will be considered in the equations to follow. The terms airplane 
moment and tail load and airplane inertia will refer to one -half of the 
airplane only. 

The theoretical aerodynamic wing loads on one semispan are given by 
equation (24) of reference 4 as 


M - q * s H ^ [ £ 


0 n 

1 


IS! + Si 


where 


i«g} 




(Bl) 


includes the wing-root angle of attack, which was not 


included in the 




of reference 4. The geometric and downwash 


matrices of this equation were calculated for a vortex system of nine 
horseshoe vortices per semispan. The values of the section lift-curve 
slopes rriQ were determined through the use of equation (Bl) as applied 
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to the wind-tunnel data of reference 10 and includes correction to the 
Mach number by the use of the Prandtl-Glauert equation. The required 
slopes of the nacelle pitching moment plotted against angle of attack 
were obtained from reference 12. 

In the determination of the wing structural twist, the wing was 
treated as a cantilever fastened at the intersection of the fuselage wall 
and the wing 38-percent-chord line. The wing -bending and torsional- 
stiffness distributions of references 7 and 8 were used in the determina- 
tion of the wing twist matrices . The influence coefficients of reference 1 
were used to determine the deflection matrices. 

For the fuselage load calculations, the fuselage was divided into 
seven segments forward and seven segments rearward of the wing root. The 
only loads considered on the fuselage in the section blanketed by the 
wing root were those induced from the wing. 

The loads induced on each section of the fuselage by one wing semi- 
span may be deduced from equation (A4) of reference 5 as 



where the elements of 



are expressed by 


(y 2 - h 2 - H 2 )a 2 

+ 2H 2 (y2 + h 2 ) + (y 2 - h 2 ) 2 


(B2) 


The point of application of the fuselage loads is at the same longitudinal 
position as the section loads on the wing. 

The loads on one -half the fuselage due to the fuselage angle of 
attack is given by equation (l) of reference 6 as 


l = q — sin a 
dx 


and for small angles of attack 


D 
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when the fuselage is considered in sections this expression becomes, for 
the part of the fuselage on one side of the fuselage center line 



r°n 

Ax 



(B3) 


where for incremental angles 0^ = 0^. 


The elements of j^Ko, j are thus section values of q The section 

values of A were obtained by completing the fuselage -section drawings of 
reference 9 by consideration of the line drawings of the same reference 
and numerically integrating for the desired areas. The effect of the nose 
radar dome and canopy was neglected but small modifications were made for 
the angle -of -attack nose cone. 

The load on each section of the fuselage due to wing wash from one 
semi span is given by 



(B 4 ) 


where 



was calculated by the methods of reference 4 , appendix E, 


and pertains only to one semispan. 


The sum of the loads at the two wing -attachment pins due to the 
wing quarter-chord airloads is expressed uy 


+ 


- L 2h j{ l v} 


(B5) 


Airplane Pitching -Moment Component Equations 

When the moment effects of the loads given by equations (Bl) through 
(B5) are combined with the moments due to the nacelles and the moments 
due to the elevation of the line of action of the wing forces in the 
x-direction as caused by the wing deflections, the final airplane moment 
equations are established. The moments due to the fuselage loads are 
obtained by inserting the appropriate moment-arm matrices into the load 
equations . 
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Fuselage pitching -moment contributions . - The airplane moment due to 
the loads induced on the fuselage by one wing semispan is given from 
equation (B2) as: 




(B6) 


The 

fuselage 


airplane moment contribution of each fuselage section due to 
angle of attack is given from equation (B3) as 



r O-j 

r ° | 

Ax 

K a 



(B7) 


where for incremental angles 


Oj, = a f 


The airplane moment contributed by each fuselage section due to wing 
wash of one semispan acting on the fuselage is given from equation (b 4) as 



rO- 

" O- 

r O -1 

Sc 

Ax 

Xx, 




(B8) 


Wing pitching -moment contributions .- The airplane moment due to the 
wing quarter-chord airload is given by 


m = 




rO-j r o 1 
2h 


'W 


(B9) 


The airplane moment due to the line of action of the chordwise 
component of wing section drag being elevated by the deflection of the 
wing is given approximately by 



(BIO) 


NACA RM L56LO5 


27 


where cos a is assumed, equal to 1 and 

f d = ( c d,o + ^) q - 2hc 

The airplane moment due to the line of action of the chordwise com- 
ponent of the lift vector being elevated by the deflection of the wing is 
given approximately by 


M- 

r° i 

r°i 

- -1 


- 0- 

hr 

Be 

s 2 

w - 

K 


where sin a is assumed equal to a 
the wing-section angle of attack. 


^jf-ojor - [hr\ < B11 ) 

and [j 3 2 jjz w j + -jhoja*. + is 


Nacelle pitching -moment contripution . - The moment due to the nacelle - 
thrust line being elevated by the deflection of the wing as caused by the 
0.25-chord airloads is given by 


{m} = Ti[k a ]|z v | + T 2 [k b J{z w } (B12) 

The moment due to the nacelle -thrust line being changed by the inertia 
of the wing is given by 

m = T^&j^ + T 2 ^i,2 (B13) 

The airplane moment due to the nacelle pitching moment caused by the 
wing twist resulting from the wing 0.25-chord airloads is given by 


m 






(B 14 ) 


The airplane moment resulting from the nacelle pitching moment due 
to wing-root angle of attack and wing inertia twist is given by 


m = 


i^ m a,e) 1 + (Si,e) 2 ( < l* Sca r) + ( C m^ e) x (^* Sca i,l) 


H*J 2 ( q * Sc_a ^) 


(B15) 
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Load Solutions 


r°i 


-0 1 

r 0 — 

_ r D 

> 

Dg L and 

_ L w_ 


In equations (BIO) and (Bll), the terms 

occur. Since it is necessary to have all of the moment equations involve 
only the unknowns < and a preliminary solution of equation (Bl) 

l > r~0 r O-i 


must be made in order to obtain approximate values for 
r o 


% 


‘D 


, and. 


Lyj. This is accomplished by solving equation (Bl) simultaneously with 


a preliminary normal-force equation and letting oj. and 




be the 


unknowns. In the preliminary normal -force equation, the loading -j.nboard 
of the fuselage wall was assumed to be that of the most inboard wing 
section. Thus, preliminary force equation in matrix form was 


|2h, 2h, 2h, . . . 2h, 2h + jl) = nw or 




Whether to use nw or Lm • depends on whether the solution is for 

x 

normal -acceleration or pitching -acceleration loads. From these approxi- 
mate wing loads, approximate wing deflections are determined by taking 

inertia effects into account. Thus, approximate values are also 

obtained . 

The moment contributions of equations (BIO) and (Bll) are small and 
could be neglected without appreciable error but were retained in the 
solutions of this paper as a matter of interest and solution refinement. 
If these equations were neglected, then the need for a preliminary solu- 
tion would be eliminated since the purpose of the preliminary solution is 
to reduce equations (BIO) and (Bll) from second and third power expres- 
sions of the wing loads to an approximate linear form. Calculations show 
that these nonlinear effects do not appreciably affect the assumptions of 
linear variations of the loads at the wing -attachment pins as assumed in 
the data analysis of this report. 

The solution of the preliminary equations provide excellent data to 
keep constant check on the final solution. 

Equations (B13) and (BI5) may be added to form the equation: 


m = eoj. + f 


(Bl 6) 
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Also after values of |f d | and JlJ and 

the remaining moment equations of (b 6) to (Bll) may be 
equation: 



are established 
added to form the 


H ■ MM + { h K + { J } 


(B 17 ) 


where each row of each matrix represents the moment contribution of a 
section of the fuselage or the wing. The elements of each column of 
equation (BI 7 ) may be added to form a moment equation 


m = l_Pj|Z w | + Bo^ + G 

to which equation (Bl6^ may be added to form a moment equation of the 
same form 


m 


L«J{*v} + + V 


Also adding the tail contribution yields the total airplane aerodynamic- 
moment equation 


m 


L«j{ ! w} 


Ra^. 


+ V + L+Z 


tH 


(B18) 


Similar additions of the force equations will yield a similar 
equation 


F = 



+ WOL^ + L.J- 


(B19) 


Solving equations (Bl), ( Bl8 ) , and (B19) simultaneously will yield the 
wing loads , the tail load, and the wing-root angle of attack. The values 

of |a,gj> and the summation of the moments and the forces must be pre- 
scribed for the flight condition for which the loads are to be determined. 
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Solution for the Fuselage Deflections Due to the 


Normal Acceleration Loads 

In the calculation of the deflection due to the normal -acceleration 
loads, the aerodynamic loads were obtained by solving simultaneously 
equations (Bl), (Bl 8 ), and (B19). The force equation was set equal to 
the n = 1 inertia load and the moment equation was set equal to zero. 

The values of jag| were given by 


^| a g j r n _-^ “ j 01 ! , n| ( ® ) + ^0 j | a rj' + ) 


The simultaneous solution gives values of Oj., L^., and l w . With 
these knowns, the fuselage loads other than and the attachment-pin 

reactions were obtained from equations (B2), (B 3 ), and (BA). The loads 
at the attachment pins were determined from the solutions of the equations 


F fs x + F rs x = m nac + m z ^ w + mi^ w+riac + m D + m L ^ c + imp 
F FS + F RS = hr + L i,w+nac 


(B21) 


The terms on the right-hand side of the equations can be obtained by 
using the previously established moment and load equations and the now 

known values of anc ^ -| a rj’* 

The fuselage inertia loads were obtained from the dead -weight dis- 
tribution. The distributions of fuselage weights, fuselage stiffness, 
and the variations of the fuel-tank centers of gravity with the amount 
of fuel are given in figures 12, 13, and 1^, respectively. The fuel-tank 
centers of gravity were obtained from the data given in reference 13 - 

From the known fuselage aerodynamic and inertia loads, the bending 
moments were obtained by numerical integration of the equation 



l dx dx 


(B22) 
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The deflections were obtained by the double numerical integration of the 
equation 


Z’ 



dx 


(B23) 


The fuselage weight and stiffness distributions were obtained from 
reference 9 and from data obtained from the manufacturer, being modi- 
fied in the manner described in the text. The stiffness distribution 
used was for the B-V7B airplane. The differences between the stiffnesses 
of the B-47A and B-47B fuselages are believed to be small. 


Solution for the Deflections Due to 
Pitching Acceleration 

In the calculation of the deflections due to pitching acceleration, 
the procedure was similar to that used for the normal-acceleration 
deflections. In this case the force equation was set equal to zero and 
the moment equation was set equal to the inertia moment. The values of 

j a g| vere defined by 

{ a g}4=l = + [ S o]{ a r} + (B24) 


The fuselage inertia loads were obtained from the dead-weight dis- 
tribution and consideration of the accelerations at the various fuselage 
sections . 
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TABLE I-. FLIGHT CONDITIONS 


Flight 

Run 

Mach 

number 

Airplane 

weight, 

lb 

Altitude, 

ft 

Dynamic 
pressure, 
lb/sq ft 

Center of 
gravity, 
percent M.A.C. 

Fuel forward 
auxiliary tank, 
lb 

Fuel forward 
main tank, 
lb 

Fuel center 
auxiliary tank, 
lb 

Fuel center 
main tank, 
lb 

Fuel aft 
main tank, 
lb 

4 

19 

0.70 

108,900 

25,000 

265 

20.97 

0 

7,700 

0 

12,500 

7,200 

4 

20 

.60 

108,700 

25,000 

189 

20.91 

0 

7,600 

0 

12,500 

7,ioo 

4 

21 

.51 

108,400 

25,000 

128 

20.85 

0 

7,500 

0 

12,400 

7,000 

6 

12 

.78 

108,700 

30,000 

267 

13.10 

1,500 

9,200 

375 

10,050 

6,100 

6 

13 

.74 

108,400 

30,000 

244 

13.12 

1,500 

9,000 

375 

10,050 

6,000 

6 

14 

.69 

108,200 

30,000 

215 

13.16 

1,500 

8,950 

375 

10,050 

6,000 

6 

15 

.64 

107,600 

30,000 

187 

13.03 

1,500 

8,700 

375 

10,050 

5,800 

10 

3 

.60 

127,200 

30,000 

159 

22.64 

1,300 

13,300 

1,300 

15,500 

14,300 

10 

5 

.68 

126,300 

30,000 

200 

22.39 

1,300 

13,100 

1,300 

15,100 

14,000 

10 

6 

.72 

126,100 

30,000 

230 

22.51 

1,300 

13,000 

1,300 

15,000 

14,000 

10 

7 

.76 

125,400 

30,000 

255 

22.96 

1,300 

12,600 

1,300 

14,700 

14,000 

10 

9 

.80 

124,900 

30,000 

275 

23.28 

1,300 

12,300 

1,300 

14,500 

14,000 
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TABLE II.- EXPERIMENTAL FUSELAGE DEFLECTION COEFFICIENTS 


(a) Deflection coefficients due to 
zero lift plus droop, Z 0 + Zj_ 


Flight 

Run 

Target 

1 

2 

3 

4 

19 

2.41 

2.82 

3.81 

4 

20 

2.17 

2.44 

3.41 

4 

21 

2.03 

2.22 

2.74 

6 

12 

2.33 

2.63 

2.98 

6 

13 

2.00 

2.10 

2.52 

6 

14 

1.94 

2.16 

2.47 

6 

15 

1.89 

1.89 

1.93 

10 

3 

2.12 

2.32 

2.21 

10 

5 

2.58 

2.76 

3.12 

10 

6 

2.59 

3.12 

3.47 

10 

7 

2.73 

3.07 

3.76 

10 

9 

2.71 

3.25 

4.25 


(c) Deflection coefficients due to 
pitching -acceleration loads, 


Flight 

Run 

Target 

1 

2 

3 

4 

19 

0.83 

1.85 

5.38 

4 

20 

1.40 

2.73 

6.14 

4 

21 

I.58 

2.52 

6.24 

6 

12 

1.14 

2.07 

4.78 

6 

13 

1.80 

2.82 

5.48 

6 

14 

1.65 

2.79 

5.34 

6 

15 

1.23 

3.22 

5.44 

10 

3 

3.12 

4.30 

9.52 

10 

5 

2.00 

3.69 

7.84 

10 

6 

1.71 

2.38 

7.40 

10 

7 

1.39 

2.69 

7.48 

10 

9 

1.54 

2.85 

6.37 


(b) Deflection coefficients due to 
normal -acceleration loads, Zn 


Flight 

Run 

Target 

1 

2 

3 

4 

19 

1.76 

2.66 

5.01 

4 

20 

1.93 

2.86 

5.07 

4 

21 

1.97 

2.93 

5.42 

6 

12 

1.82 

2.76 

5-33 

6 

13 

2.00 

3.06 

5.40 

6 

14 

1.97 

2.88 

5.19 

6 

15 

1.93 

3.01 

5-44 

10 

3 

2.64 

3.77 

6.62 

10 

5 

2.24 

3.43 

5.96 1 

10 

6 

2.25 

3.10 

5.82 

10 

7 

2.15 

3.27 

5.77 

10 

9 

2.24 

3.26 

5.58 


(d) Deflection coefficients due to 
pitching-velocity loads, Z^ 


Flight 

Run 

Target 

1 

2 

3 

4 

19 

-0.17 

- 0.84 

-2.11 

4 

20 

-.62 

-1.26 

-.67 

4 

21 

-.20 

-.47 

-33 

6 

12 

-1.16 

-.48 

•72 

6 

13 

-1.37 

- 2.12 

-I.71 

6 

14 

-.16 

-•95 

-.92 

6 

15 

-1.17 

-1.28 

-1.82 

10 

3 

-1.50 

-1.96 

- 3-41 

10 

5 

-.82 

- 1.32 

-2.72 

10 

6 

-.99 

-.97 

- 3.00 

10 

7 

-.83 

-1.81 

-2.75 

10 

9 

-1.60 

- 1.94 

- 3.94 
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TABLE III.- TYPICAL TARGET DEFLECTIONS AND STANDARD ERRORS 

[Flight 10, Run 7] 


Target 

Target-deflection coefficient and standard error for - 

s z 

( Z ° + Z i) 


z 4 

Z d 

1 

±0.027 

2.735 ±0.037 

2.153 ± O.O57 

1.595 t 0.169 

-0.852 ± 0.249 

2 

±0.022 

3.077 ± O.OJO 

5.272 ± 0.050 

2.696 ± 0.158 

-1.818 ± 0.205 

3 

to .056 

3.763 ± 0.076 

5.778 ± O.O77 

7.482 ± 0.551 

-2.756 ± 0.518 
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TABLE IV.- COMPARISON OF SOME QUANTITIES OBTAINED IN THEORETICAL FUSELAGE 
CALCULATIONS WITH EXPERIMENTAL RESULTS FROM SEVERAL SOURCES 


Condition 

Method 

Weight, 

lb 

Angle of 
attack, 
da/d n 
(a) 

Tail load, 
dL t /dn 

(b) 

Fuselage station 
of center of 
gravity 

Wing root 
shear at 
y = 54 
(c) 

Pitch inertia 
Iy, lb -in. -sec 

(1) 

Flight 

4 

Run 

20 

Experimental 

108,700 

O.O85 

902 

(+L t up) 

617.4 

46,000 

15,881,000 

Calculated 

109,072 

.084 

1,372 

613.4 

49,000 

1 ^, 659,749 

Flight 

10 

Run 

7 

Experimental 

125,400 

.064 

2,231 

620.6 

52,300 

18,320,000 

Calculated 

125,772 

.067 

2,916 

617.2 

55,800 

17 , 734,298 


(^Experimental data from reference 5. 

( ^Experimental data from preliminary results of strain-gage measurements on the wing and tail 
of the test, airplane as obtained in the Langley Flight Research Division. 

( c ) 

Experimental data from preliminary results of strain-gage measurements on the wing and tail 
of the test airplane as obtained in the Langley Flight Research Division. 
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Figure 1.- Photograph, of test airplane. L-86692 
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Figure 2.- Three views of test airplane. 
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; at fuselage sta. 1*60 
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sta. 1*02 
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sta. 835 


Figure 3.- Top and side view of test airplane showing optigraph camera and fuel-tank locations. 
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Figure 4.- Topical time history of fuselage deflections. 
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Lines are deflections due to zero-lift + droop 
+ normal load factor as determined from analysis 
of the data . 
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.- Fuselage deflections measured at various load factors in a 

typical inn. 
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Figure 6.- Typical fuselage-deflection curves obtained from the deflection 
coefficients, q* = 265 lb/sq, ft; weight = 108,900 lb; M = 0.70. 
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Figure 7*- Variations of the fuselage-deflection coefficients for the 

various runs. 
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Figure 8.- Comparison of measured fuselage deflection with that calculated 

using the deflection coefficients. 
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Figure 9.- Theoretical fuselage bending moments and comparison between theoretical and experimental 
deflections due to the normal-acceleration loads. Weight = 125,^00 lb; M = O. 76 ; q.* = 255 lb/sq ft 
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Figure 10.- Theoretical fuselage bending moments and comparison between theoretical and experimental 
deflection due to the pitching-acceleration loads. M = O.76; q* = 255 lb/sq ft; weight = 125,4-00 lb. 
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Figure 11.- Comparison of experimental and theoretical fuselage deflections due to the normal- 

acceleration and the pitching-acceleration loads. M = 0.60j weight = 108,700 Ibj 1* = 189 lb/sq ft. 
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Figure 12.- Approximated weight distribution of the fuselage including weight of the pilots, 

instruments, and unavailable fuel. 
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Figure 13.- Fuselage bending-stiffness distribution for positive bending (downward). 
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Figure 14 .- Variation of fuel tank center of gravity with amount of fuel (data from ref. 12). 
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